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ABSTRACT 

We present a new method to use 24 /Km observations to estimate the total infrared luminosity, 
L(TIR), and the star formation rate of star- forming galaxies across the redshift range < z < 3. This 
is accomplished by parameterizing the spectral energy distribution (SED) of a galaxy as a function of 
the L(TIR) surface density, I]^(xir)5 instead of the traditional L(TIR). Our reformulation allows an 
accurate description of star-forming galaxy SEDs, including the aromatic emissions, out to redshift 
z r-^ 3.0. We test the new 24 fim L(TIR) indicator against L(TIR) measured from stacked far-IR 
photometry at redshift < z < 3. We show that a monochromatic 24 /j,m observation can be used 
to estimate L(TIR) consistent with the values determined from the multi-band far-IR measurements 
such as those from Herschel, on average within 0.1— dex. This approach allows use of observed-frame 
24 /im observation to determine the SFR of star-forming galaxies accurately across the currently 
expected peak of the star formation history of the universe. The success of this method indicates that 
a majority of IR-luminous star-forming galaxies in the universe are not strongly nuclear concentrated 
as are the local merger-induced LIRGs and ULIRGs. 

Subject headings: cosmology: observation — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies 



1. INTRODUCTION 

Mid-Infrared (Mid-IR) wavelengths provide a unique 
window to study the evolution of galaxies across the 
history of the universe, particularly at redshift 1 < 
z < 3 where the s tar formation rate fSFR) of the uni- 
verse peaks (e.g., iHopkins fc Beacoin[l2006[ ). At red- 
shift z > 0.5 to 3, a majority of star formation took 
place in obscured environments and most of the UV 
photons from hot young stars, the signpost for star for- 
mati on, are rep rocessed by dust into IR emission (see 
e.g. iLe Floc'h et al. 2005t iPerez-Gonzalez et all [20051 : 
iDole et all l2006t iBuat et all I2007D . The opticallv ex- 
tincted environment, especially in galaxies with IR lu- 
minosity, L(TIR), greater than the characteristic lumi- 
nosity, L* , poses inherent challenges for optical and UV 
determinations of the SFR. 

There are two approaches to study the SFR using IR 
observations. Both involve using IR photometry as a 
tracer of L(TIR), which in turn traces the number of 
UV photons from short-lived massive stars and allows the 
SFR to be inferred (see e.g. Kennicutt 1998). The first 
approach is to use multi-band IR photometry covering 
the wavelengths of the peak of the IR emission to con- 
strain the L(TIR) directly. Recent examples use the Her- 
schel Space Observat ory to observe star-forming galaxie s 
at 100 - 500 Atm (e.g. lElbaz et al.|[20Tol: iRex et al.i[20lol ). 
The complete characterization of the spectral peak of the 
dust emission provides a good measurement of L(TIR) 
and SFR. However, the sample size of this approach is 
compromised by requiring detection of an object in mul- 
tiple bands that, towards longer wavel engths, are increa s- 
ingly limited by confusion noise fe.g. lDole et al.|[2004[) . 
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The second approach is to use a monochromatic IR lu - 
minosity to trace L(TIR) and SFR. iCalzetti et al.l ()2007D 
have shown that the rest-frame 24 /um emission is well 
correlated with the ionizing photons that directly trace 
star formation, while shorter wavelengths such as 8 /im 
have a larger scatter due to the emission from stochas- 
tically heated aromatic molecules. Monochromatic far- 
IR indicators (70 /im and 160 /im) are affected by the 
emission from colder dust heated by old stars, which re- 
sults in an increasing scatter in the SFR calibration as 
a function o f wavelength, particul arly in low-luminosity 
IR galaxies. ICalzetti et al.l (|2010a[) found that the corre- 
lation between 70 and 160 /xm emission and SFR suffers 
from a scatter larger than t hat at 24 /im (jRieke et al.l 
120091: iKennicutt et al.l I2OO90 . These factors render 24 
/im a particularly suitable monochromatic SFR indica- 
tor locally. 

Beyond the local universe, however, the observed 24 
/im band probes wavelengths in the spectral energy dis- 
tribution (SED) of a galaxy that are increasingly af- 
fected by aromatic emissions (e.g. polycyclic aromatic 
hydrocarbons or PAH). At z ~ 2 the observed 24 //m 
band probes the strongest ar omatic emission co mplexes 
at 6.2, 7.7, and 8.6 /im fe.g. ISmith et al.l 120071 ). which 
helps boost the 24 /im flux and thus aids detection of the 
galaxy at the cost of greater uncertainty in the estima- 
tion of the SFR. However, the most challenging issue at 
higher redshifts lies not in the intrinsic scatter, but rather 
that the aromatic emission in high-redshift star-forming 
galaxies at a given L(TIR) appears to be significantly 
strong er relative to £(TIR) than in their local counter- 
parts ()Rigbv et al.ll2008HFarrah et al.|[200llTakagi et all 
l2010f ). A direct consequence of the apparent evolution of 
aromatic strength is that a mid-IR SFR indicator will 
overestimate i(TIR) and hence SFR by up to a factor 
of 5 — 7 if local SED templates are used for k-correction 
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without reconsidering the aromatic strength relative to 
the template L(TIR). This is reported as the "mid-IR 
excess" problem in recen t Herschel results ("Elba z et alj 
l2010HNordon et alJ[20Tol ). The evolution of the aromatic 
features in star-forming galaxy SEDs must be taken into 
account in using the mid-IR to indicate L(TIR) and SFR. 

A motivation to use 24 /xm observations to study the 
SFR despite the availability of the direct and accurate 
L(TIR) and SFR measurements from Herschel is that 24 
fiia observation still detect more individual galaxies and 
probe weaker SFR at any given r edshift from < z < 3 
(see Figure 4 of E lbaz et al.li2011| ) . In addition to the the 
Spitzer 24 /im data in very deep survey fields, the large 
areas surveyed by Spitzer (e.g., SWIRE and the Bootes 
field) will not be fully surveyed by Herschel to the same 
depth in terms of L(TIR) or SFR. The wavelength is 
also uniquely accessible with current and future space 
missions such as WISE (2009, 3-25 ^.xn) and JWST 
(-2017, 0.6 - 28 /im). 

Although the intrinsic scatter of the aromatic emission 
remains an issue, the larger systematic effect of the red- 
shift evolution of aromatic emissions can be addressed by 
building on our previous result (Rujopakarn et al. 2011) 
that the L(TIR) surface density of a star-forming galaxy, 
^L(TiR): serves as a good indicator of the SED. This pa- 
per demonstrates how this result can help select the ap- 
propriate SED template for star-forming galaxies beyond 
the local universe. We derive a new L(TIR) and SFR es- 
timator using the observed 24 /xm fi ux over the redshift 
range of < z < 3 and based on the iRieke et"al1 ()2009f ) 
SED library for local galaxies. 

The paper is organized as follows. Section [2] discusses 
the evolution of the SEDs of star-forming galaxies and 
the use of El(tir) as a guide to select the appropriate 
SED for high-redshift galaxies. The derivation of the new 
monochromatic L(TIR) indicator is described in Section 
13.11 We compare the new L(TIR) indicator with other 
independent measurements as a test in Section 13.21 and 
present the final L(TIR) and SFR indicator with further 
tests using larger samples in Section |4l Section [5] dis- 
cusses the validity of rest-frame 8 — 24 /im as a tracer 
of L(TIR) and SFR, along with some further implica- 
tions of our approach. Throughout this paper, we as- 
sume a ACDM cosmology with fi^ = 0.3, f^A = 0.7, and 
Hq = 70 km s~"'^Mpc~^. We will refer to galaxies with 
L(TIR) in the range of 10" -10^^ and those with L(TIR) 
> 10^^ i© as Luminous Infrared Galaxies (LIRGs) and 
Ultraluminous Infrared Galaxies (ULIRGs) , respectively, 
or collectively U /L IRGs. We followed the convention of 
Rieke et al.l (I2009D a nd ad opt the definition of L(TIR) of 



Section [3] 

2.1. Sl(tir) o,s an Indicator of SED for Star-Forming 

Galaxies 

Until now, most estimations of the L(TIR) and SFR 
for moderate and high redshift (z > 0.5) galaxies using 
monochromatic IR photometry rely on an assumption 
that the SED of a star-forming galaxy does not evolve 
with redshift. SED templates of IR galaxies are con- 
structed from local galaxies, which is an inhomogeneous 
population comprised of both normal disk star-forming 
galaxies and those with nuclear star-formation that ap- 
pears to have been induced by galaxy interactions. The 
most IR-luminous galaxies found locally invariably rep- 
resent int eraction-induced sta r formation to various de- 
grees (e.g. iVeilleux et al"1l2002t ). The consequence is that 
when we apply local IR-luminous SED templates to de- 
termine the SFR of high-redshift galaxies we make an 
implicit assumption that the high-redshift star-forming 
galaxies are similarly dominated by interaction-induced 
starbursts. The validity of this assumption becomes crit- 
ical at z > 0.5 where the IR luminosity density of the uni- 
verse is dominated by U/LIRGs in the luminosity range 
at which, locally, IR galaxies are virtually all dominated 
by sub-kpc nuclear regions in interacting galaxy systems. 

Observations have shown that actively star-forming 
galaxies at high redshift are different from their local 
counterparts in at least three major aspects. First, 
the IR SEDs at high-redshift exhibit colder dust tem- 

f ieratures, T^^, than local galaxie s at fixed L(TIR) 
Pope et all 120061: iSvmeonidis et all 120091: [ivRTzzin et al.l 
|2010[) . In other words, these galaxies have dust temper- 
atures consistent with local galaxies with lower L(TIR). 
Second is that the sizes of the star-forming regions in 
IR galaxies at high-redshift are substa ntially larger than 
those of local galaxies at fixed i(TIR). iRuiopakarn et afl 
(|201lD found that IR-emitting galaxies of all luminosi- 
ties at intermediate and high redshift have similar phys- 
ical IR-emitting size, mostly on the order of 4 — 10 
kpc, similar to local normal moderate-luminosity star- 
forming galaxies. By contrast, the local U/LIRGs, 
from which the SED templates for their correspond- 
ing luminosity are constructed, are significantly smaller 
with diameters on the sub-kpc scale (see Figure 3 of 
IRuiopakarn et al.l 120 lit ). The extended size of luminous 
star-forming galaxies at high redshift implies a structure 
similar to that of local normal star-forming galaxies and 
differing from local U/LIRGs despite their luminosity in 
the U/LIRG range. This common structure is recog- 
nized in recent literature as the " main sequence" of star- 
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with L(TIR) obtained by integrating the SED template 
over 5 — 1000 /xm. 

2. QUANTIFYING THE IR SED EVOLUTION 

In this section, we will discuss the challenges of using 
mid-IR observations, specifically Spitzer MIPS 24 /im 
photometry, to estimate L(TIR) and SFR both locally 
and at high redshift. Relevant issues are the determina- 
tion of an overall approach for using Sl(tir) as an in- 
dicator of a galaxy SED, and the choice of SED library. 
We will then apply this approach to describe appropri- 
ately the SED of galaxies at high redshift in terms of the 
local SED templates, which will be discussed in detail in 
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2010t lElbaz et al.ll20ial20lH ). where star-formation oc 
curs at a relatively steady rate rather than in bursts, and 
in disks or clumps, rather than in merger nuclei. In its 
simplest form, this picture suggests that these galaxies 
differ among themselves (and from local lower luminosity 
star-forming galaxies) primarily in the SFR surface den- 
sity, SsFR, and hence in IR luminosity surface density, 

^ L(TIR)- 

iChakrabarti fc McKed (|2005[ ) modeled the radiative 
transfer in centrally heated dusty sources and found a 
general relation among the luminosity-to-mass ratio, the 
surface density, and the shape of the SED. They pre- 
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Fig. 1. — The relationships of the L(TIR) surface density, 
^L(TlR)i S'li'i iJ(TIR) differ for local U/LIRGs (circles), from which 
we construct the SED templates for IR galaxies, and the U/LIRGs 
beyond the local universe (grey squares). Upward triangles rep- 
resent lower limits for both the local and high-redshift galaxies 
whose S^z-piR) is derived from CO size measurements that could 
be systematically larger t han the rest of the samp le based on radio 
continuum measurement lIRuiopak arn et al. 20110 . The dotted line 
and the dashed line are the fits for high-redshift and local galaxies 
(ignoring CO lower limits), respectively. Ruiop akarn et al.[ II2011I) 
show that galaxy SED characteristics are indicated by Si(TiR)- 
This allows one to describe high-redshift galaxy SEDs with appro- 
priate local SED templates. For example, a high-redshift galaxy 
with L(TIR) of lO^'''^ LqwIU have aromatic emissions and absorp- 
tion features consistent with a local galaxy with the same Sl(tir) > 
indicated by an arrow, in this case a local galaxy with L(TIR) of 
~ 10^^'^ Lq. The method to "map" local SED te mpla tes to high- 
redshift galaxies is described in detail in Section 13.11 The solid 
line is a local fit yielding the final L(TIR) calibrate d to the FIR 
observations from BLAST and Herschel (see Section [XT} . 



diet that the behavior of SED templates for loeal galax- 
ies might not extend to high redshift. Extended sizes 
and colder SEDs were shown to be connected throug h 
a Td-^L{TiR) relation explored by iChanial et al.l (|2007l ). 
who theoretically described the IR-emitting region of 
a star-forming galaxy as an isothermal cloud optically 
thick to optica l wave lengths and optically thin in the IR. 
IChanial et all (|2007f ) also show that L(TIR), T^, and 
the extent of the IR-emitting region are related by a fun- 
damental plane, of which the L{TIR)-Td relation (e.g. 
ISoifer et al.lll987l ) is one manifestation. 

Third is the aforementioned evolution of the strength 
of the aromatic features, which grow stronger at high 
redshifts at a fixed L(TIR). Quantitatively, for instance, 
iRigbv et al.l (pOOl ) found a galaxy with L(TIR) of lO"-^ 
Lq at z = 2.5 to exhibit aromatic emission features 
consistent with a local galaxy with an order of mag- 
nitude lower L(TIR). The appropriate description of 
aromatic feature strength evolution is crucial because 
the overestimation of the resulting L(TIR) would other- 
wise exceed even the uncertainties due to intrinsic scat- 
ter (~ 0.3— dex; iCalzetti et a l. 20071 especially where 
the aromatic con i plexes redshift into 24 /im band. In 
iRuiopakarn et al.l ()20lH ). we have demonstrated that the 
underlying parameter explaining the apparent evolution 
is the Ei(TiR). 

In Figure m adapted from Figure 4 of IRuiopakarn et aLl 



(|201lD . we show Sl(tir) as a function of L(TIR) along 
with trend lines describing the relationship in each pop- 
ulation, one for local galaxies with L(TIR) above 10^"'^ 
Lq and another for high-redshift galaxies. The surface 
area of the IR-emitting region used to estimate Sl(tir) is 
measured by radio continuum (e.g., 1.4 GHz), Paschen- 
a, or 24 /Ltm observations such that they reflect the size of 
the star-forming regions but are insensitive to the contri- 
bution from old stellar populations. The main sequence 
can be drawn in this figure from the local normal star- 
forming galaxies (non-U/LIRGs) onto the high-redshift 
star-forming galaxies. Both local and high-redshift trend 
lines agree at low luminosity, but the trend that in- 
cludes local U/LIRGs diverges from the main sequence at 
i(TIR) ^ 10^^ Lq where the local galaxy sample starts 
to be drawn from the interaction-dominated U/LIRGs 
populations. Drawing a horizontal line at fixed Sl(tir) 
in Figure |T] for the IRigbv et al.l ()2008D example, we ex- 
pect the galaxy on the main sequence with L(TIR) of 
10^^-^ Lq to have spectral features consistent with a 
galaxy on the local sequence with L(TIR) of ~ 10^°-^ 
Lq, which agrees with t he observed aromatic features. 
IRuiopakarn et al.l (|201lD further demonstrated the use 
of Sl(tir) to predict SED behavior by predicting the 
24 /im-to-1.4 GHz flux ratios and matching the observed 
stacked aromatic spectrum of high-redshift ULIRGs to a 
local SED template. 

We therefore adopt the assumption that Si(TiR,) is 
the dominant parameter controlling the average SED. 
Given the large range of ^^(xir) represented in both the 
main sequence and the local U/LIRGs, it is then possible 
to choose an appropriate S ED to represen t high- redshift 
U/LIRGs. Apart from the IChanial et all (|2007| ) funda- 
mental relationship, this approach is also physically mo- 
tivated in that the aromatic emission emerges from the 
outer surfaces of the photodissociation regions (PDRs; 
e.g. lTielensll2008l ). so the surface density of the star for- 
mation rate in a galaxy is among the controlling factors 
for the aromatic-emitting surface of a galaxy. 



2.2. Selection of the Local Reference for the Galaxy SED 

The most widely used lo cal SED librar i es are those 
of ICharv & Elbazl (l200l) . IDale fc Heloul (l2Q0l . and 
iRieke et al.l ( 2009r ) . The former two libraries were based 
on spectra taken by ISO that do not fully span the 6 — 24 
/im wavelength range that the 24 fim band probes at red- 
shift z = 3 — 0. The recent infrared spectrograph (IRS) 
observati ons from Spitze r, cove ring 5.2 — 38 /xm, show 
that the I Charv fc Elbai (|200lD SED templates (here- 
after, CEOl) have suppressed aromatic features at very 
high luminosity and do not sufficiently take into ac- 
count the strong silicate absorption features at '^ 10 
/im, resulting in weaker aromatic bands but stronger 
overall emissi on in the mid-IR inconsistent with actual 
galaxy SEDs (lArmus et al.ll2007l : lRleke et al.ll2009t) . The 
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IDale fc Heloul (|2002l' ) SED library is optimized to de- 
scribe moderately-luminous local star-forming galaxies 
with i(TIR) < 10^^ Lq and hence the lack of silicate 
absorption becomes more evident, especially at higher 
L(TIR) with larg er extinction. We therefore adopt the 
IRieke et al.l (|2009l) SED library (hereafter, R09) for a lo- 
cal SED reference. 
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iRieke et al.l ()2009D constructed SED templates sepa- 
rately for normal star-forming galaxies (sub-LIRG) and 
those in U/LIRG range in order for the final SED li- 
brary to cover L(TIR) of lO^-^^ to lO^-OO Lq. The 
U/LIRGs SED templates were constructed from a sam- 
ple of 11 local LIRGs and ULIRGs with high quality 
Spitzer IRS spectra as well as photometric data covering 
optical to r a dio w avelengths (see Figure 1, 2, and 3 of 
IRieke et alJ ()2009[ )1 These U/LIRGs were chosen such 
that their IR emission is dominated by star-forming ac- 
tivity. Construction of the SED library is done in two 
steps. First, these 11 galaxies were used as a basis to 
assemble 11 archetypal SED templates spanning 0.4 /Ltm 
to 30 cm wavelength. Their IRS spectra (5 — 38 fim) 
were joined to the photospheric emission and far-IR dust 
emission components in a series of tests to ensure both 
spectral continuity and appropriate flux calibration. Sec- 
ond, these archetypal templates are combined with dif- 
ferent weights to produce the final averaged SED tem- 
plates. The template weights are optimized by match- 
ing synthesized IR colors from the combined template 
to IR colors of observed galaxies with known L(TIR) 
from the IRAS REGS (Sanders et al. 2003 ). The R09 li- 
brary construction was extended to star-fo rming galaxies 
with s ub-LIRG L(TIR) by combining the IDale fc Heloul 
(|2002D SED library with the recent mid-IR spec t ral li- 
brary based on IRS observations from lSmith et al.l (|2007f ) 
using the same IR color fitting technique as the U/LIRGs 
template construction. The use of IR color to help guide 
the combination of archetypal templates (i.e. fitting 25/8 
/im, 25/12 //m, and 60/25 fim colors simultaneously) 
helps ensure that the final templates represent the av- 
erage property of real galaxies even though they are con- 
structed from a limited sample. The R09 SED library 
construction is de scribed in detail in the Appendix of 
IRieke et al.l (|2009l ). 

3. A MID-IR ESTIMATOR FOR L{TIR) AT < Z < 3 

The construction of a 24 ^m SFR indicator has two 
steps: (1) to construct an i(TIR) estimator using 24 fim 
photometry, and then (2) to derive a relationship be- 
tween L(TIR) and SFR. We discuss construction of the 
L(TIR) estimator in Section [3Tl and test it by compar- 
ing our L(TIR) estimates to i(TIR) derived from other 
independent measurements in Section 13.21 The test is 
then used to refine our L(TIR) indicator, which we pro- 
vide in final form in Section [J] along with the relationship 
between L(TIR) and SFR. 



3.1. Construction of a Mid-IR L(TIR) Indicator 

The process of calculating the i(TIR) typically in- 
volves finding the best-match SED for a set of ob- 
served flux es at a given redshift and then integrating the 
SED (e.g. 'Perez-Gonzalez eLalj 120051: iLe Floc'h et~all 
120051 : Magnclli et al. 2009). For monochromatic indi- 
cators, this means fitting an SED using a single data 
point, whic h leav es significant room for uncertainties. 
IRieke et"all (|2009| ) circumvent this process by first as- 
signing a value of L(TIR) to each SED template and 
then determining their corresponding SFR as well as 
the monochromatic luminosity through the desired fil- 
ter. For the 24 fiui band, in particular, the assignment 
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formal fit to the L(TIR)-E^(xiR) relationship of local galaxies (Sec- 
tion [3TTJ; Col 3: resulting L(TIR)ncw associated to SED template 
shape from the formal fit; Col 4: stretch factors from the final 
L(TIR)-S£(XIR) parameterization calibrated to FIR observations 
(Section 14.11 1: Col 5: resulting L(TIR)new from the final parame- 
terization. 

t The IRieke et al.l ||2009| ') SED templates with L(TIR)oid > lO^^ 
Lq are omitted because their L(TIR)new are higher than lumi- 
nosities of observed galaxies. These SED template shapes are thus 
unlikely to represent real galaxies at high redshift. 

of the 24 /im luminosity of a template follows the ratio 
of L(24 ^m)-to-i(TIR) based on the local IRAS data 
(see Fig. 8, Fig. 15, and formula A6 of IRieke et all 
|2009| ). Once each template has an associated 24 /xm lu- 
minosity, it is possible to calculate the k-corrections and 
subsequently the expected monochromatic flux for each 
template as a function of redshift. For L^, oc vf,y, the k- 
correction, Kcorr{z), and the flux, f^{24:fim), are related 

by 



Kc 



,{z) = (l + z) 



fi,{v == (1 + 2:)j/obs) 



/^(24^m) 



A-kDl /i/,obs 



Lu,rcst{'24:fim) 



i^24 



-tvcorr \Z j 



(1) 



(2) 



where D^ is the luminosity distance for an object at 
redshift z. The 47rDL^/24,obs for this equation has the 
unit of Jy cm^. The relationship between the observed 
monochromatic flux from each template and the tem- 
plate's monochromatic luminosity at a given redshift is 
approximately linear, which allows for a linear fit to be 
performed at each redshift (i.e. for a redshift grid) to 
determine a set of coefficients that convert the observed 
flux at a given redshift to the luminosity or any desired 
quantity associated with the template (i.e. L(TIR) or 
SFR). For example, a relationship between L(TIR) and 
the observed 24 /im flux, /24,obs7 would have the form 

logL(TIR)^i^ = Aiz) + B{z) [log(47ri?iV24,obs) - C] , 

(3) 
where A{z) and B{z) are the intercept and the slope from 
the linear fit, respectively, C is a zero-point to reduce 
covariance in the fit parameters, and L(TIR)oid refers to 
the L(TIR) associ ated to each of the R09 SED templates. 
IRieke et"al] (|2009[ ) use this method to tabulate the coef- 
ficients to convert the monochromatic fiuxes in various 
bands to the SFR of galaxies. 
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Because of the evolution of the IR SEDs, a differ- 
ent approach is needed to derive s u itable templates for 
them. Following 'R uiopakarn et ahl ()2011[ ). an appropri- 
ate choice of SED to use for k-correction beyond the local 
universe is the one corresponding to the same Sl(tir) lo- 
cally. This is equivalent to increasing the L(TIR) associ- 
ated with each of the R09 SED templates by a factor de- 
termined by the ratio of luminosities with equal Sl(tir) 
on the local U/LIRGs trend to those on the main se- 
quence. Wc will refer to this ratio as the stretching factor, 
Si, for each SED template. Note from Figure[l]that Si is 
nearly negligible at X(TIR) ^ 10^^ Lq and reaches three 
orders of magnitude at L(TIR) of 10^* Lq. To deter- 
mine Si quantitatively, we parameterize the relationship 
by fitting a 2'^'^-order polynomial to the local U/LIRG 
relation (excluding the lower limits of Sl(tir) from CO 
observations) and a linear fit to the main sequence, and 
then take the ratio of i(TIR) on the main sequence fit 
to the R09 template L(TIR) at the same Sl(tir) (see 
Figure [TJ the Si values are tabulated in Tabled]). Since 
we effectively increase the luminosity of each template 
by a factor of Si, the observed flux will also be increased 
by the same factor. Equation |3] can thus be rewritten as 



logL(TIR)„„^ = log[SiL(TIR)„i,] 
A'{z) + B'{z) [\og{AnDL'S,f2A^o\.s) ~ C] 



(4) 



The new set of coefficients, A'{z), B'{z), and C" that 
use the appropriate SED templates for main sequence 
star-forming galaxies at redshift < z < 3 can be de- 
termined by re-fitting equation |4l We have limited the 
fitting range to only encompass the "stretched" luminos- 
ity of logL(TIR)ncw < 14. This is because the Si are 
very large at high luminosity and the stretched luminosi- 
ties are well above the luminosity range occupied by real 
galaxies (the 10^'* Lq cut off is chosen because we expect 
this to be at or beyond the upper luminosity limit of star- 
forming galaxies). For a given template, the shape as a 
function of redshift remains the same. The effect of the 
Si is to change the spacing between templates. The re- 
lationship between luminosity and observed flux remains 
well approximated as linear, with residuals < 0.05 dex. 

The fitting coefficients, A'{z) and B'{z), to relate the 
observed Spitzer 24 /im flux to i(TIR) as a function of 
redshift, are shown in Figure [2] (as dotted lines). The 
zero-point, C", for the formula is 45. Thus the relation 
is 



logL(TIR)„^„ = A'{z) + B'{z) [\og{A7TDL'f24, 



obsj - 45] , 

(5) 

As is the case for lRieke et al.l ()2009D . the feature in A'{z) 

and B'{z) at z '^ 1.4 is due to silicate absorption, which 

depends on template luminosity. 

3.2. Comparison of L{TIR) to Far-IR i(TIR) 

Measurements 

We test our L(TIR) indicator in two steps. In this 
section, we apply it to small samples of galaxies with 
accurately measured i(TIR) from far-IR bands covering 
the spectral peak of dust emission (e.g. 200 — 500 /xm). 
Using the result from this step, we make fine adjustments 
to the indicator in Section 14.11 and perform a series of 
blind tests on large samples with more limited far-IR 




Redshift (z) 

Fig. 2. — The coefficients of the fits for a relationship between the 
observed Spitzer 24 /im flux and the total IR luminosity (Equation 
O as a function of redshift. The coefficients A'{z) and B'{z) are 
shown by the dotted red and blue lines. The red and blue solid 
lines represent the final coefficients A* (z) and B* (z) resulting from 
a fit calibrated to the FIR L(TIR) measurements. The latter set 
of coefficients are tabulated in Table [2] 

photometric coverage (e.g. 24 — 160 fim) to investigate 
the performance of the indicator. 

For the flrst step, we compared the L(TIR) estimated 
from our method with those derived from the integration 
of SEDs fltted to the FIR observations from Herschel as 
well as the Balloon-borne Large Aperture Submillimeter 
Telescope (BLAST; e. g. iPascale et ahl I2Q08I: iDve et al.l 
[20091: lEaleset al.ll2009t) . We also compared our L(TIR) 
to those estimated using the CEOl and the original R09 
SED libraries. 

The Herschel data used for this t est were taken from 
the Herschel Lensing Survey (HLS; lEgami et al.l 120101 1. 
The sample consists of 19 galaxies located behind the 
Bulle t Cluster at redshifts 0.4 < z < 3.2 ()Rex et al.l 
l2010f) . These sources are detected in at least two Her- 
scfte/ bands (at 100—500 fj,m) and many are also observed 
in LABOCA 870 m icron and AzTEC 1.1 nim ma ps of the 
fleld (Wilson et~all[200l [Johansson et al.ll2010 t). These 
measurements tightly constrain the peak of the FIR SED 
and therefore provide accurate estimates of L(TIR). We 
have excluded one object from this sample, IILS18, be- 
cause its large lensing magniflcati on (54 x) is not well 
constrained due to nearby objects (jRex et al.ll2010( ). 

The BLAST data we used arc a subset of 27 galax- 
ies from the photometric catalog of BLAST sources lo- 
cated in the Ch andra Deep Field South presented in 
iDve et ahl ("2009") that have spectroscopic redshift mea- 
surements (Ealcs et al. 2009). The BLAST photometry 
was originally calibrated bas ed on o bservations of the 
red supergiant star VY CMa (|Truch et al...20091 Recent 
SPIRE FTS observations of this calibration source reveal 
significant emission lines within the BLAST bandpasses 
(jRover et al.ll201Cil ). Convolving the observed spectrum 
for VY CMa with the BLAST fllter transmission curves 
indicates excess flux at the level of 0.06— dex compared 
with the originally assumed blackbody. For this work the 
BLAST fluxes have therefore been rescaled by this fac- 
tor to accurately match the SPIRE measurements. The 
galaxies in this sample have photometric measurements 
in at least two of the BLAST bands (at 250 - 500 /xm) 
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which are augmented in most cases by Spitzer 24, 70 and 
160 micron data for SED fitting. The BLAST data are 
ultimately limited by source confusion at submillimeter 
wavelengths; the catalog is therefore dominated by main 
sequence galaxies at redshift z < 0.5. All of the higher 
redshift galaxies are contaminated b y AGN and are ex - 
cluded based on their classification bv lEales et al.l ()2009D . 

Figure [3] compares L(TIR) estimated from this work 
with that derived from the HLS and BLAST data, as well 
as the same comparison for i(TIR) based on the CEOl 
and R09 SED libraries. The figure demonstrates the 
"mid-IR excess" by which the 24 /im monochromatic in- 
dicator over-estimates L(TIR) compared with that mea- 
sured by integrating the SED fitted to far-IR data. The 
mid-IR excess is significant for the monochromatic 24 
/im-derived _L(TIR) using both the CEOl and the original 
R09 SED libraries. However, the excess is alleviated by 
using our adjusted indicator and the resulting monochro- 
matic 24 ^m-derived L(TIR) agrees well with the far-IR 
observations. 

It appears from Figure [3] that the CEOl SED library 
is less prone to overestimating L(TIR) for high- redshift 
galaxies compared to the R09 library. A careful compari- 
son of the two template sets has, however, shown that the 
omission of silicate absorption features in the CEOl SED 
templates results in an overestimation of aromatic-region 
flux in the local SED templates that coincidentally helps 
partially compensate the PAH-region flux evolution at 
high redshift. The R09 templates correctly include the 
silicate absorption locally, and hence suffer a larger over- 
estimation of L(TI R) at intermediate and high redshift 
(see the appendix of lRieke et al.l (|2009[ ) for a detailed dis- 
cussion of uncertainties in high-redshift behavior of the 
SED templates) . We have now shown that the R09 tem- 
plates can be used to extrapolate for L(TIR) accurately 
from monochromatic 24 fim observations given an appro- 
priate choice of SED template based on the Sl(tir)- 

4. FINAL APPROACH TO DETERMINING L(TIR) AND SFR 

In this section we discuss a slightly modified calibration 
of our template selection to provide the final indicator 
along with the relationship between L(TIR) and SFR. 
The final indicator is then tested against large samples of 
star-forming galaxies from cosmological surveys to verify 
its integrity over < z < 3. 

4.1. An Alternate Calibration Based on a Modified 
Parameterization to the Local L(TIR)-E/^(xiR) 
Relationship 

Figure |3] hints at a small (< 0.2— dex) yet systematic of 
underestimation of L(TIR) at high luminosity, especially 
at L(TIR) > 10^^ Lq. a possible cause of the underes- 
timation is that there are other underlying higher-order 
parameters that affect the SED shape at high luminosity 
other than our simple assumption that the SED shape 
is indicated by I]^(xir). This necessitates future explo- 
rations of the SED dependence on a broader range of 
parameters, such as redshift, luminosity, metallicity, and 
environment. A second possibility leading to the under- 
estimation could be the uncertainties in the local fit of 
the L(TIR)-E^(xiR) relationship that is constrained by 
a sample of 44 local galaxies, only 25 of which are in 
the U/LIRG luminosity range. A larger sample of lo- 
cal galaxies in the U/LIRGs range would be ideal, but 
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Fig. 3. — Co mparison of L(TIR) der ived from this work and the 
SED libraries of lCharv fc Elbaj pOOH') and lRieke eFall ||2009|'| . la- 
beled as CEOl and R09, respectively, to the Z/(TIR) measured by 
integrating the far-IR multi-band photometry from Herschel Space 
Observatory and the Balloon-borne Large Aperture Submil lime ter 
Telescope (BLAST). The mid-IR excess problem ^Section I3.2|l is 
shown to be alleviated by the new L(TIR) indicator. Our 24 ^in- 
based L(TIR) is derived from the local L(TIR)-Sspp relation fit 
that excludes CO lower-limit measurements (Fig. [TJ. Note that in 
the bottom plot, the averaged data point is shown at the median 
luminosity in the corresponding bin. At high-L(TIR), this indica- 
tor exhibits a slight underestimation of L(TIR) compared to the 
FIR measurements. Nevertheless, our monochromatic estimate of 
L(TIR) is consistent with the FIR L(TIR) measurements within 
0.2-dcx on average. This underestimation could be a result of the 
uncertainties in the parameterization of the local L(TIR)-Sspr 
relationship as well as the possibility of higher-order underlying 
parameters other than the S^z-piR) that could affect the shape 
of high-redshift galaxy SEDs. We provide a separate relationship 
based on the same ^i^itir) rnethod that improves agreement of 
our monochromatic L(TIR) estimates and the FIR L(TIR) mea- 
surements as shown in Fig. |4] 

the sample size is limited by the number of known local 
U/LIRGs that do not harbor significant AGN activity 
and have good size measurements. 

In light of the uncertainties of the local L(TIR)- 
^L(TiR) relationship (dashed line in Figure [T]), it is pos- 
sible to modify the parameterization of this relationship 
by calibrating the monochromatic L(TIR) to the FIR 
L(TIR) observations from Herschel and BLAST. This 
is done by optimizing a 3"^-order polynomial describing 
the local L(TIR)-I]sfr relationship to minimize the dif- 
ferences between the resulting L(TIR) estimated by our 
indicator and the i(TIR) based on FIR measurements. 
From this parameterization, we calculate new stretching 
factors, S* , tabulated in Table [I] which results in the 
final coefficients for equation [51 A*{z) and B*{z), tab- 
ulated in Table [2] and shown in Figure [2] in comparison 
with A'{z) and B'{z) derived in Section |3T] The modi- 
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Fig. 4. — L(TIR) estimated from 24 fim observations using the 
^L(TIR) ^s S'li indicator of SED shape and a modification of the 
local L{TIR)-Sjj(xiR) relation to calibrate the final L(TIR) to the 
FIR L(TIR) measurements. The symbols are the same as in Fig. 
\3\ This relationship allows one to estimate L(TIR) and SFR con- 
sistent with the FIR measurements (e.g. Herschel and BLAST) on 
average within 0.1— dex using monochromatic 24 /ivci observations. 

fied relationship (solid line in Figure[T]) is a less steep rela- 
tion that bends toward the CO lower-limit measurements 
of Sl(tir)- The final _L(TIR) estimates agree with the 
FIR -L(TIR) measurements on average within 0.1— dex 
(Figure HI). 

With the final parameterization, we also tabulated the 
luminosity of the appropriate R09 SED template that 
should be used to describe star-forming galaxies given 
the 24 fiva observed flux and redshift in Table [31 Even 
at the bright-end of the flux range at high redshift (e.g. 
/24 — 3.0 mJy at z = 3.0, which corresponds to L(TIR) 
of 10^^-^ i©), the appropriate R09 template for such 
galaxy is that of a local ULIRG with i(TIR) of IQi^'' 
Lq. In fact, it is evident from the table that most IR- 
luminous star-forming galaxies at < z < 3 exhibit spec- 
tral characteristics of local galaxies with L(TIR) in the 
range of 10" -5 x 10" Lq. 



4.2. Testing the Final L(TIR) Indicator with Samples 
from Large Cosmological Surveys 

There are two aspects of the test required to establish 
the reliability of our L(TIR) indicator in real-life usage, 
based on estimating L(TIR) in blind samples of galaxies 
from cosmological surveys. First is to quantify the scat- 
ter of the resulting L(TIR) compared to L(TIR) derived 





TABLE 


2 


Coefficients of 


THE 


Final Fits for 


THE Relation 


Between Spitzer 24 


flM 


Flux and the 




L{TIR 




z 


A*(z) 


B*{z) 


0.0 


2.669 


0.973 


0.2 


2.303 


1.026 


0.4 


1.960 


1.070 


0.6 


1.915 


1.080 


0.8 


1.490 


1.118 


1.0 


1.014 


1.168 


1.2 


0.775 


1.209 


1.4 


0.675 


1.231 


1.6 


0.725 


1.206 


1.8 


0.669 


1.185 


2.0 


0.633 


1.184 


2.2 


0.597 


1.189 


2.4 


0.527 


1.200 


2.6 


0.429 


1.222 


2.8 


0.315 


1.245 


3.0 


0.232 


1.263 



Note. — A*{z) 
and B*{z) are the co- 
efficients for equation 
[5] derived from the fi- 
nal parameterization 
of the relationship 
calibrated to FIR ob- 
servations (solid line 
in Fig. 1). See de tails 
in Section [3.1l and r4.1l 

from multi-band SED fitting. Second, and more impor- 
tantly, is to verify that the indicator accurately describes 
the SED features of galaxies over < z < 3, i.e. that the 
mid-IR excess problem has been resolved. We conclude 
from the tests in this section that our L(TIR) indicator 
accurately estimates L(TIR) within 0.1— dex on average 
compared to L(TIR) determined from SED- fitting and 
with scatter of less than 0.26— dex. 

We adopt three samples for these tests: (1) the 
70 /xm-selected COSM OS sample (n = 751) from 
iKartaltepe et al.l (|2010[ ): (2) the stacked 70 /z m and 160 
Atm ph otometry of 24 /im-selected sample from lLee et al.l 
(I2OIOI) in COSMOS [n > 35,000); and (3) the stacked 
24 /xm, 70 ^m, 160 ^m, 1.4 GHz, and 610 MHz obser- 
vations of a NIR-selected sample {n = 3, 172) from the 
Extended Chandra Deep Field-South survey (ECDFS) 
(jBourneet all 120111) . Unlike the Herschel and BLAST 
samples, these L(TIR) estimates rely on extrapolating a 
fitted-SED to the dust emission peak. The goal of this 
section is thus not to re-calibrate our indicator, but to 
test it on a large sample across a broad redshift range. 

To quantify the scatter, we compare our final 
L(TIR) calibration to the i(TIR) estimates from 
the 7 /im-selected COSMOS sample (jKartaltepe et al.l 
I2OIOI ): the entire sam ple contains 150 3 galax ies at 
0.0 < z < 3.5. IKartaltepe et al.l ()2010f ) esti- 
mated L(TIR) by fitting IR photometry from Spitzer 
at the 8, 24, 70, and 160 /im bands to an ex- 
tensive collection of SED libraries a nd selected the 
best fi t via y^ minim i zation among the iCharv fc Elba3 
(•MTl-'Dale fc Helou' 1*200^ JLagache et al.l (|2003D. and 
iSiebcnmorgen fc Kriigel (,20()7l ) templates. For our com- 
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cal L(TIR)— S^criR) relationship lj4.ip is used here. 

parison, we exclude sources with X-ray luminosity > 
10'*^ erg/s, those with radio-excess, and those exhibit- 
ing power-law SEDs in the Spitzer IRAC bands to 
minimize AGN contam ination. Above 10^^'^ Lq, the 
IKartaltepe et al.l (|201C0 sample contain s very few sources 
that do not harbor AGN or QSO. IKartaltepe et alj 
(|201ClD show that sources with 160 /xm detections are 
biased toward objects with colder dust, whose L(TIR) 
may not directly correlate with SFR, and therefore we 
also exclude these sources. Lastly, we limit our com- 
parison sample to those with uncertainties in SED-fitted 
L(TIR) < 0.35 dex to avoid comparing to outliers with 
large uncertainties in luminosity. After applying these 
cuts, we have 751 sources left. These sources have red- 
shifts ranging from z = 0.07 to 1.81 with a mean and 
median redshift of 0.52 and 0.43. The luminosity of this 
subsample ranges from L(TIR) of 10^-^ to 10^^-^ Lq and 
both mean and median L(TIR) of lO^-^ and 10"-^ Lq, 
respectively. 

The comparison shown in Figure [5] indicates an offset 
of ;^ 0.05— dex betw e en ou r L(TIR) estimates and those 
of IKartaltepe et all ()2010[ ) across the luminosity range. 
This offset is iden tica l for L(TIR) es t hxiate s based on 
IRieke et all (l200l and lCharv fc Elbazi ([2002) at L(TIR) 
< 10"'^"'^ Lq and is likely due to a difference in convention 
of L(TIR) between the indicators. Apart from this offset, 
the estimates are consistent within 0.1— dex at L(TIR) 
of 10^'^ — 10^^'^ Lq with an average internal scatter of 
0.26— dex. Some fraction of this scatter must ari se in the 
assignment of _L(TIR) bv lKartaltepe et "aLJ (|2010[ ): hence, 
it is an upper limit relative to our L(TIR) estimates. In 
ad dition to estimating the scatter, the comparison with 
the IKartaltepe et al.l ()2010D sample allows us to test our 
L(TIR) estimates particularly at z ^ 0.4 where the 17 
fixn aromatic complex redshi fts int o the 2 4 /im band. The 
agreement with Kartaltcpe et al.l (|2010[ ) indicates that 
our SED assignment method adequately describes the 
aromatic emission strength. 

Apart from testing the 24 ^m monochromatic L(TIR) 
indicator, we also use the modified local L(TIR)-Ssfr 
relationship to derive the A{z) and B{z) coefficients to 
estimate L(TIR) from the 70 ^m observed flux (e.g. 
Spitzer MIPS or Herschel PACS). The resulting coef- 
ficients yield monochromatic 70 fim i(TIR) consistent 
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Fig. 5. — The scatter of our L(TIR) indi cator is estimated by 
comparing our final calibration (Section 14. 1 [I to those from multi - 
band SED fitting in the COSMOS sample lKartaltepe et"aI|[20Tol) 
of 751 galaxies at redshift 0.0 < z < 1.8. Red circles, diamonds, 
and grey circles represent L(TIR) estimates from this work, the 
R09, and the CEOl SED library, respectively. Color-coded points 
are individual galaxies. Error bars refiect the internal scatter in 
each luminosity bin, which is on average 0.26— dex. Apart from the 
systematic offset of ~ 0.05— dex across the luminosity range, our 
monochromatic L(TIR) indicator yields results consistent within 
0.1-dex of multi-band SED fitting. The agreement of L{TIR) at 
redshift z ~ 0.4 and 2; ~ 1 where the 17 ^m and 11 fim aromatic 
emission complexes redshift into the 24 fim band suggests that our 
method accurately predicts the spectral features. 



with the values from the original IRieke et"al1 ()2009f ) 
estimator a t L(T IR) < 10^^ Loand overestimate the 
IRieke et al.l ()2009[ ) L(TIR) by < 0.15-dex at higher lu- 
minosities. We do not anticipate strong evolution of the 
SED in the rest-frame wavelength range probed by the 
observed 70 /xm band at redshift z = — 3 (70 — 17.5 
fim) because it is still in the Wien tail region of the dust 
blackbody emission, unlike the 24 /im band. Thus the 
Rieke et al. (2009) coefficients can be used to estimate 

L(TIR) and SFR from 70 fiui flux. 

The second test employs the iLee et all (|2010[ ) sample 
that stacks 70 /zm and 160 fj,m. observations of over 35,000 
galaxies selected at 24 ^m in COSMOS. The stacks were 
done in bins of 24 /xni flux (0.06 < /24 < 3.00 mJy) 
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and redshift (0 < z < 3). We note that although 
iKartaltepe at alj (|2010[ ) found that individual sources de- 
tected at 160 /xm are biased toward galaxies with cold 
dust as mentioned earlier, the same study also found that 
the stacked 160 /im flux from galaxies not detected indi- 
vidually is not affected by this bias and hence the 160 /im 
stacked photometry can be used reliably. Lee et al. then 
fit SEDs to the average 24 /^m, 70 ^m and 160 /^m fluxes 
in each stacked b i n wit h the same procedure as that of 
IKartaltepe et alJ ()2010( ) to estimate L(TIR), which we 
use to test our indicator. We have excluded bins where 
the fraction of X-ray sources (probabl y associat e d wit h 
AGN) exceeds 0.1 based on Figure 3 of llee et all (|2010D . 
For the actual L(TIR) calculation we use the average 
value of 24 /xm flux and redshift in each bin (N. Lee, 
2011, private communication) to calculate the average 
L(TIR) of the bin, whi c h can be compared with the es- 
timates from Lee et al.l ()201Qil . The result from this test 
is presented in Figure [HI 

Fina lly we used the stacking result from lBourne et all 
(|2011l ). Their sample of 3172 galaxies was selected using 
Spitzer ITiAC 3.6 /im and 4.5 /im photometry (i.e. prob- 
ing stellar mass) in ECDFS. The selection of the sam- 
ple in the near-IR provides an independent sample that 
might reveal selection effects (if any) inherent to mid- 
and far-IR selections (e.g. the Lee et al. and Kartaltepe 
et al. samples). Bourne et al. then stacked observations 
of the sample in 7 redshift bins from < z < 2.0. These 
observations include 24 /im, 70 /im, 160 /im, 1.4 GHz and 
610 MHz; the stacked fluxes of the Spitzer MIPS bands 
were then used to fit a M51 SED to estimate L(TIR), 
which we use to compare with our indica to r (Fig ure Et. 

The comparison with iLee et al] (|2010D and 
iBourne et al.l (|2011l ) shows excellent agreement and 
indicates that our description of galaxy SEDs appropri- 
ately describes the observations out to z ^ 3. There is 
no indication of systematic trends emerging at z '^ 1 and 
z ^ 2 where the 11 /im and 7 /im aromatic complexes 
redshift into the 24 /im band. 

4.3. L{TIR)-SFR Relation 

Finally, SFRs can be determined by making use of 
the relationship between L(TIR) and the rest- frame 
L(24 /im) and subsequentl y the L(24 /im) and SFR, 
given by iRieke et al.l (|2OO90 with small modifications. 
The introduction of the stretching factor, which is ef- 
fectively to re-normalize the SED templates, requires 
a modification of the relationship between L(TIR) and 
L(2A iMn). The o riginal fit as given in equation A6 of 
IRieke et all (12001 ) is 



logL(TIR)oid = 1.445 + 0.945 logL(24) 



(6) 



The new relationship can be derived by re-fitting equa- 
tion [6] with the stretching factor, 5^, multiplying both 
L(TIR) and L(24 /im) for each template i. The re-fitted 
relation is given by 



logL(TIR), 



1.113 + 0.980 logL(24) (7) 



To determ ine the SFR froin L(24 /im), the calibra- 
tion given bv lRieke et al.l ()2009D remains valid. However, 
that calibration has a term that corrects for a decrease 
in L(24 /im)/L(TIR) above L(TIR) = 10" Lq. Since 
the correction is motivated by an increase of the optical 
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Fig. 6. — Comparison of L(TIR) from our monochromatic 24 
fim indicator to L(TIR) measured by stacking far-IR observations. 
The comparison samples are from Lee et al. (2010) who determined 
L(TIR) at < 2: < 3 by stacking 70 fim and 160 ^m obser v ations 
of over 35,000 COSMOS galaxies, and from B ourne et al.l ll20Tll) 
whose L(TIR) was determined from stacking a broad range of ob- 
servation from 24 fim to 610 MHz of 3,172 galaxies at < z < 2 
in the ECDFS-FIDEL sa mple. Open a nd filled symbols represent 
L(TIR) estimates from the lRieke eFalj (2009) SED library and this 
work, respectively. The agreement between estimated and mea- 
sured L(TIR) also implies that a majority of star-forming galaxies 
at these redshift resides in the "main sequence" of star-forming 
galaxies. 

depth at high L(TIR) that prevents the mid-IR emis- 
sion from escaping, the threshold at which optical depth 
becomes significant depends directly on the geometry of 
the galaxy. In the same way that the extended struc- 
ture of the galaxy affects the IR-emitting environment, 
the optical depth will consequently be lower for a given 
L(TIR) and the luminosity threshold where the optical 
depth should apply has to be scaled up by a stretch- 
ing factor as well. The Si corresponding to the original 
threshold is a factor of 3.3, yielding a luminosity thresh- 
old of 3.3 X 10^^ Lq. This is equivalent to L(24 /im) of 
4.1 X 10^° Lq. Therefore the relationship between SFR 
and L(24 //m) is given by 

SFR(Mq yr^i) = 7.8 x 10"^" L(24/im, L©) 

for 5x 10^ Lq < L(TIR) < 3.3x 10" Lq or 6x 10* Lq 
< L(24 /im) < 4.1 x lO^o Lq. For L(TIR) > 3.3 x 10" 
Lq or L(24 fim) > 4.1 x lO^o Lr. 



^0j 



SFR(M0 yr^i) = 7.8 x 10"^" L(24/iTO, L©) 



[2.41 X 10"" L(24/im, Lq)] 



0.048 



(8) 
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This calibr ation is based on that of iKennicutti (|1998[ ) 
but with the iKroupal ()2002D IMF, which yields S FRs a 
factor of 0.66 of those assuming the lSalpeteii (|1955[ ) IMF. 
A library of IDL routines for the L(TIR) and SFR indi- 
cator is available at our websitcQ. These relations apply 
at z> 0.3; within the local volume t he L(TIR) and SFR 
relations given by iRieke etlUI ()2009[ ) should be used. 

5. DISCUSSION 

Here we discuss the validity of using 24 /im observation 
at < z < 3 to derive L(TIR) and SFR. At these red- 
shifts, this is equivalent to providing justifications that 
the rest-frame 6 — 24 /im luminosity that encompasses 
aromatic emissions can be used to estimate i(TIR). We 
will also discuss the implications of our L(TIR) indicator 
in light of the test results in Section 14.21 



5.1. Validity of the Aromatic Emission as a L(TIR) 
and SFR Indicator 

Our indicator utilizes the luminosity from observations 
that correspond to rest-frame wavelengths of 24 /im to 6 
/im at z = to 3. This necessitates an investigation 
into the validity of using the mid-IR-dominating aro- 
matic emission as a tracer for L(TIR). We discuss two 
major factors affecting the aromatic feature strength: the 
metallicity dependence and the presence of AGN. 

The correlati on of the aromatic l umin os ity to L(TIR) 
was studied by 'Rigopoulou et alj (|1999( ). IRoussel et al.l 
(|2001) a nd Elbaz et al. (2002) using ISO] later by 
Wu et al.l (I2005D with Spitzer, and recentl y by 
Elbaz et al.l (|2011[) with Herschel. iCalzetti et all |2007l ) 



used the SINGS sample to show that the scatter and the 
non-linearity of the relationship between the extinction- 
corrected 8 /im luminosity and the extinction-corrected 
Paa luminosity (which is a good tracer for SFR) is un- 
satisfactorily large for the former to serve as a SFR indi- 
cator. However, their Figure 10 suggests that the largest 
source of scatter is due to the metallicity differences 
among the subsamples. For high-metallicity systems 
(Z > V^ Zp whic h is equivalent to 12+log(0/H) > 
8.2). lCalzettil (|2010b[ ) reports that the stellar-continuum- 
subtracted aromatic emission shows a good correlation 
with the SFR. F urthermore, lEngelbracht et al.l (|2008D 
and iSmith et al.l (|2007) found that the ratio of aro- 
matic luminosity to L(TIR) does not vary significantly 
a t metallicity^, 1/3 Zcr, . 

ISmith et ahl C2007^ found that AGN can significantly 
supp r ess the aromatic emission (see also, Moorwood 
[l98a IRocheet all Il99ll : IGenzel et al.l [199811 . The sup- 
pression of aromatic luminosity is shown in Figure |7| 
(le ft) based on the tota l aromatic luminosities measured 
by ISmith et al .l (l2007f) and a metallicity measurement 
fromTMousta kas et al.i (|2010f ). AGN, however, do not 
affect the aromatic emis s ion o f the entire host galaxy. 
iDiamond-Stanic fc Riekd (|2010[ ) use Spitzer IRS to com- 
pare nuclear and non-nuclear spectra of nearby Seyfert 
galaxies and found that the aromatic features in the nu- 
clear spectra are clearly suppressed while the aromatic 
emission is of normal strength in the outer disk. There- 
fore, with regard to using aromatic emission to trace 
L(TIR): (1) the aromatic suppression effect due to AGN 

^ http://ircairLera.as.arizona.edu/rujopakarn2011 



is likely limited to the nuclear region and aromatic emis- 
sion arising from the rest of the galaxy should still pro- 
vide a good tracer of the SFR; and (2) the longer wave- 
length emission such as the 11.3 /im aromatic emission 
complex, which the 24 /tm observation probes at z ~ 1, 
could still serve as a SF R indicator despite the presenc e 
of low-luminosity AGN piamond-Stanic fc Riekell2010() . 

Galaxies at z ~ 2 have on average 0.3— d ex lower 
metallicity than local galaxies (|Erb et al.l 120061 ). There- 
fore, galaxies more massive than 3 x 10^ M© should 
be sufficiently metal rich (> 1/3 Zq) that the 
L(TIR) /aromatic calibration is valid. If we consider the 
ratio of aromatic-luminosity-to-L(TIR) as a function of 
metallicity in the left panel of Figure |7| within the range 
of metallicity expected at high redshift, the scatter of the 
relationship between aromatic luminosity and i(TIR) is 
in fact about 0.1-dex in the non-Seyfert/LINER sample. 
This suggests aromatic luminosity to be a good indicator 
for L(TIR) at high redshift given that an effort is made 
to exclude AGN from the sample (such as the AGN ex- 
clusion criteria employed in tests in Section |4.2|) . 

We further test the validity of aromatic emission (e.g. 
8 /im rest- frame observation) as a measure of L(TIR) 
in a sample of Hll-dominated galaxies with moderate 
to high metallicity in the right panel of Figure [71 In 
this Figure, the 8 /im observations are from Spitzer 
and the L(TIR) estimates are based on the IRAS all- 
sky survey using 6, 12, 25, and 60 /im observations 
([Sanders et al.ll2003[ ). This sample is chosen such that 
it mimics the population of star-forming galaxies ex- 
pected at high redshift: high-metallicity normal star- 
forming galaxies without AGN. The scatter of the re- 
lationship between L(8 /im) and i(TIR) in Figure [7| is 
0.14-dex w i th th e ratio L(TIR)/L(8 /im) of 4.3 ± 1.6. 
lElbaz et al.l ([201 1[ ) report a similar ratio measurement lo- 
cally (L(TIR)/L(8 /xm) = 4.8^1:7) and also show that the 
ratio is universal for high-redshift galaxies in GOODS- 
Herschel, further reinforcing the applicability of aromatic 
emission as an L(TIR) and SFR indicator. 

5.2. The Eddington Luminosity of ULIRGs 

It has been suggested that local ULIRGs are opti- 
cally thick to mid-IR photons and radiation pressure 
may play a role in lim i ting their maximum luminosity 
(iTho mpson et al.ll2005[ : [Younger eralll2008t [ThompsonI 
2009). If the degree of compactness is assumed to be 
similar for both the local and high-redshift ULIRGs, an 
obvious problem would arise that the IR emission of high- 
redshift ULIRGs with L(TIR) commonly found far above 
their local counterparts would exceed the Eddington lu- 
minosity. Our results ease the concerns about the Ed- 
dington limit of ULIRGs at high redshift for two rea- 
sons. First, galaxies at higher L(TIR) are affected by 
the mid-IR excess to a greater extent. Our new indi- 
cator would therefore lower the L(TIR) estimates espe- 
cially for high-luminosity galaxies. To test the effects of 
L(TIR) reduction to the luminosity functions of galax- 
ies, we construct a set of IR LFs of 24 /im sources at 
< 2: < 2.5 using the 24 /xm observations from COS- 
MOS (e.g. [Sanders et al[l2007l ) and find that the bright- 
end of the LFs in all redshift bins are below L(TIR) of 
~ 10^^ L© (quantitatively, galaxies with L(TIR) above 
lO^'^ Lq are rarer than 10~^/Mpc^/dex in all redshift 
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Fig. 7. — (Left) The sources of scatter in using aromatic emissions (e.g. PAH) to predict L(TIR) and SFR are primarily the metallicity 
dependence and the presence of AGN. Blue dots and red stars show the ratios of aromatic luminosity to L(TIR) as a function of metallicity for 
HII and AGN dominated local galaxies, respectively (Smith ct al. 2007). For Hll-dominated galaxies with metallicity greater than 1/3 Zq 
(solid vertical line), which is similar to the enviromnent of high-redshift star-forming galaxies, the correlation of aromatic luminosity, LpAHi 
and L(TIR) has a ~ 0.1— dex scatter, suggesting that aromatic emission could serve as a good L(TIR) and SFR indicators at high redshift. 
(Right) The correlation between the rest-frame L(8 ^m) and L(TIR) for local star-forming galaxies that are dominated by Hll-emission 
and have high metallicity, similar to those expected at high redshift, has a scatter of 0.14— dex. The small scatter underlines the potential 
of using aromatic emission as a L(TIR) and SFR indicator. 



bins), which suggests that the maximum L(TIR) of star- 
forming galaxies is in general less than 10"'^'^ Lq. Second, 
the succes s of the corrections for t he size evolution de- 
scribed in iRuiopakarn et al.l (|201lD as well as this work 
to account for the general IR SED evolution indicate that 
the great majority of high-redshift IR galaxies are phys- 
ically extended, which lowers the optical depth of the 
galaxy. Both the lower L(TIR) and extended structure 
imply that high-redshift ULIRGs are emitting below the 
Eddington limit and that the maximum L(TIR) of these 
galaxies is not governed by radiation pressure. 

5.3. Modes of Star Formation 

In deriving the stretching factors from the local 
Sl(tir)"^(TIR) relation to the main sequence relation, 
there is an implicit assumption that there are two distinct 
relations: one for local galaxies with L(TIR) > 10^""^ Lq 
generally from compact nuclear regions (starburst mode) 
and another for the main sequence galaxies. This is mo- 
tivated by the fact that local U/LIRGs are rare objects 
found only by large- volume surveys such as the ffl^ 5 all- 
sky survey, and that most of them are interacting systems 
with compact central starbursts. Beyond the local uni- 
verse, other surveys do not encompass enough volume to 
find similar systems in large numbers. 

At intermediate and high redshift, the precise fraction 
of galaxies that form stars in starburst mode remains a 
subject of controversy. A tracer for starburst mode star- 
formation is the galaxy merger fracti on, which c an be 
estimated from morphological studies. iLotz et aD ()2008l ) 
report a relatively constant major merger (i.e., mass ra- 
tio > 1 : 3) fraction of 10% ± 2% for all galaxies in 
their sample at 0.2 < z < 1.2 and ~ 15% for U/LIRGs. 
iShi et ai] ()2009() find that the fraction of galaxies exhibit- 
ing asymmetric morphologies in the same redshift range 
is 20% - 3 % an d as many as ~ 50% for LIRGs. The 
ILotz et all (|2008[ ) figures represent a conservative esti- 
mate by l imitin g mergers to major mergers while the 
iShi et al] (|2009[ ) provides an upper limit, since not all 
asymmetric galaxies a re mergers. For ULIRG and LIRG 
populations at z < 1, iKartaltepe et al.| (|2010[ ) find that 
the merger fractions are ~ 50% - 80% and ~ 25% - 40%, 



respectively. However, morphological studies are increas- 
ingly compromised by the limited angular resolution and 
surface brightness sensitivity at higher redshifts. More 
importantly, disturbed morphologies found at high red- 
shift are not necessary a result of mergers but could 
also be due to clumpy or asymmetric star-formation in 
ga s-rich environments suc h as those recently observed 
bvlSwinbank et all (|2010b[ ) and iHailev-Dunsheath et all 
(|2010( ) as well as those emerging in cos mologica l simu- 
lations by lAgertz et all (2009) and Dave et al.l (|2010) . 
The range of 15% — 50% can be taken as a broad con- 
straint on the fraction of systems undergoing interac- 
tion up to z ^ 1.2, depending on how merger is de- 
fined. These merger fraction estimates imply that the 
interaction-induced starburst mode is a significant yet 
unlikely a dominant mode of star formation at interme- 
diate and high redshift. 

On the other hand, many recent observations support 
the main sequence picture of star-forming galaxies at 
intermediate and high redshift. The supporting evidence 
includes the finding that high-redshift star-forming 
galaxies lack evidence for interactions and that these 
galaxies have extended star-forming regioi is, contradict- 
ing th e picture of compact starbursts. iTargett et alJ 
()2010D find that 95% of optical morphologies of 
radio/sub-mm galaxies at z ~ 2 are consistent with 
normal disk galaxies, well described by axi-s ymmetric 
disk-like models. Similarly, I Yuma et al.l (|2011[ ) find that 
the UV morphology of BzK galaxies at z ~ 2 that have 
a single component also shows consistent Sersic indices 
with those of normal disks. Extended physical sizes 
of high-redshift star-forming galaxies, U/LIRGs, and 
submillimeter galaxies ha ve been reporte d extensively in 






the literature (iChapman et al. 2004: Muxlqw et al.ll2005l: 



Biggs fc IvisonI 120081: jBothwell et all 



i 



20091: llono e t al.'^2 009t iLehnert et all , , 

20ld iSwin bank et al.l l2010al: iTacconi et al.l 120101 



201C : ICasev et al. 



20091: iCariUiet al. 



Younger et al.ll201(]|) . In this study, the assumption that 
star-forming galaxies lie on the main sequence allows 
assignment of their SED template based on Sl(tir) 
that successfully predicts high redshift galaxy SEDs, 
and subsequently, their i(TIR). From various tests of 
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the L(TIR) indicator from this work against L(TIR) 
estimated from far-IR SED fitting for a broad range of 
galaxies out to z '-^ 3 in Section |43 the good agreement 
and smaU scatter support the initial assumption that 
a majority of star-forming galaxies beyond the local 
universe reside in the main sequence of star-forming 
galaxies. 

6. CONCLUSIONS 

The use of monochromatic mid-IR indicators (such 
as 24 /im photometry) has previously been affected by 
the differences between the structure of local U/LIRGs, 
upon which the SED templates were constructed, and 
the galaxies in the cosmological surveys being studied. 
Starbursts in the local U/LIRGs are almost invariably 
interaction-induced and very compact, while U/LIRGs 
at high-redshift are much more extended. High-redshift 
U/LIRGs typically have surface areas ^ 100 x larger than 
their local counterparts. The resulting larger surface 
area of the photodissociation regions harboring aromatic 
emission leads to stronger aromatic features at high- 
redshift mimicking an evolution in aromatic strength 
as a function of redshift. The ultimate implication of 
this morphological difference is that monochromatic mid- 
IR observations, which at 2 > 1 probe the aromatic- 
dominated SED region, will overpredict L(TIR) and 
SFR compared to far-IR measurements if the relation- 
ship of local SED templates to L(TIR) is assumed for 
k-correction. ^^^^^^^^^^^^ ^_^ 

Following the IRuiopakarn et al.l (|2011[ ) result that 
^L(TiR) can serve as a good predictor for the appropri- 
ate SED of star-forming galaxies, we determine a new 
24 /j,m monochromatic L(TIR) and SFR indicator valid 
for star-forming g alaxies at redshift 0.3 < z < 3 (the 
IRieke et"an (|2009t ) indicator should be used at z < 0.3). 



Two relationships are derived. A preliminary result is 
based on the formal fit to the L(TIR)-Ssfr relation of 
the local galaxies. We refine this fit with a small ad- 
justment to improve the agreement to the FIR measure- 
ments. The fact that the use of Sl(tir) eis an indicator 
of SED shape allows us to derive L(TIR) highly con- 
sistently wit h the FIR measurement s (within 0.2— dex) 
confirms thelRuiop akarn et al.l (|201lD result that Sl(tir) 
is indeed the dominant factor affecting the SED shapes 
of star-forming galaxies. However, the need for a small 
modification of the calibration suggests Sl(tir) may not 
be the sole parameter affecting the SED, or that there 
is a room for improvement in the sample of local star 
formation-dominated U/LIRGs to better constrain the 
local L(TIR)-S/^(TiR) relation. 

For estimating L(TIR) consistent with the FIR L(TIR) 
measurements (from, e.g., Herschel) within 0.1— dex us- 
ing monochromatic 24 //m observations, the use of the 
latter calibration (i.e. apply equation [5] with the A*{z) 
and B*{z) coefficients provided in Table [5]) is encour- 
aged. Given the L(TIR), the L(24 ^m) and SFR can 
then be determined by equations [7] and [SJ respectively. A 
monochromatic 24 /im indicator will facilitate the study 
of the SFR of galaxies at high redshifts, particularly at 
z ^ 2 where the strong 7.7 /im aromatic emission com- 
plex redshifts into the Spitzer 24 /im band as well as with 
the future 21 and 25.5 /im bands of the Mid-Infrared In- 
strument on JWST. 
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the support from the Thai Government Scholarship. 



REFERENCES 



Agertz, O., Teyssier, R., & Moore, B. 2009, MNRAS, 397, L64 

Armus, L., et al. 2007, ApJ, 656, 148 

Biggs, A. D., &: Ivison, R. J. 2008, MNRAS, 385, 893 

Bothwell, M. S., et al. 2010, MNRAS, 525 

Bourne, N., Dunne, L., Ivison, R. J., Maddox, S. J., Dickinson, 

M., & Frayer, D. T. 2011, MNRAS, 410, 1155 
Buat, v., et al. 2007, ApJS, 173, 404 
Calzetti, D., et al. 2007, ApJ, 666, 870 
Calzetti, D., et al. 2010, ApJ, 714, 1256 
Calzetti, D. 2010, arXiv: 1010.4996 
Carilli, C. L., et al. 2010, ApJ, 714, 1407 
Casey, C. M., et al. 2009, MNRAS, 399, 121 
Chakrabarti, S., & McKee, C. F. 2005, ApJ, 631, 792 
Chanial, P., Flores, H., Guiderdoni, B., Elbaz, D., Hammer, F., & 

Vigroux, L. 2007, A&A, 462, 81 
Chapman, S. C, Small, I., Windhorst, R., Muxlow, T., & Ivison, 

R. J. 2004, ApJ, 611, 732 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 
Dale, D. A. & Helou, G. 2002, ApJ, 576, 159 
Daddi, E., et al. 2010, ApJ, 714, 118 
Dave, R., Finlator, K., Oppenheimer, B. D., Fardal, M., Katz, N., 

Keres, D., & Weinberg, D. H. 2010, MNRAS, 404, 1355 
Diamond-Stanic, A. M., & Rieke, G. H. 2010, ApJ, 724, 140 
Dole, H., et al. 2004, ApJS, 154, 93 
Dole, H., et al. 2006, A&A, 451, 417 
Dye, S., et al. 2009, ApJ, 703, 285 
Eales, S., et al. 2009, ApJ, 707, 1779 
Egami, E., et al. 2010, A&A, 518, L12 
Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, H., Franceschini, 

A., Fadda, D., & Chary, R. R. 2002, A&A, 384, 848 
Elbaz, D., et al. 2010, A&A, 518, L29 
Elbaz, D., et al. 2011, arXiv:1105.253j3 



Engelbracht, C. W., Rieke, G. H., Gordon, K. D., Smith, 

J.-D. T., Werner, M. W., Moustakas, J., Willmer, C. N. A., & 

Vanzi, L. 2008, ApJ, 678, 804 
Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C, Reddy, 

N. A., & Adelberger, K. L. 2006, ApJ, 644, 813 
Farrah, D., et al. 2008, ApJ, 677, 957 
Genzel, R., et al. 1998, ApJ, 498, 579 
Genzel, R., et al. 2010, MNRAS, 407, 2091 
Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 
Hailey-Dunsheath, S., Nikola, T., Stacey, G. J., Oberst, T. E., 

Parshley, S. C, Benford, D. J., Staguhn, J. G., & Tucker, C. E. 

2010, ApJ, 714, L162 
lono, D., et al. 2009, ApJ, 695, 1537 
Johansson, D., et al. 2010, A&A, 514, A77 
Kartaltepe, J. S., et al. 2010, ApJ, 709, 572 
Kartaltepe, J. S., et al. 2010, ApJ, 721, 98 
Kennicutt, R. C. 1998, ARA&A, 36, 189 
Kennicutt, R. C, et al. 2009, ApJ, 703, 1672 
Kroupa, P. 2002, Science, 295, 82 

Lagache, G., Dole, H., & Puget, J.-L. 2003, MNRAS, 338, 555 
Lee, N., et al. 2010, ApJ, 717, 175 
Le Floc'h, E., et al. 2005, ApJ, 632, 169 
Lehnert, M. D., Nesvadba, N. P. H., Tiran, L. L., Matteo, P. D., 

van Driel, W., Douglas, L. S., Chemin, L., & Bournaud, F. 

2009, ApJ, 699, 1660 
Lotz, J. M., et al. 2008, ApJ, 672, 177 
Magnelli, B., Elbaz, D., Chary, R. R., Dickinson, M., Le Borgne, 

D., Frayer, D. T., & Willmer, C. N. A. 2009, A&A, 496, 57 
Moorwood, A. F. M. 1986, A&A, 166, 4 
Moustakas, J., Kennicutt, R. C, Jr., Tremonti, C. A., Dale, 

D. A., Smith, J.-D. T., & Calzetti, D. 2010, ApJS, 190, 233 
Muxlow, T. W. B., et al. 2005, MNRAS, 358, 1159 



Mid-IR Star Formation Rate Indicator 



13 



Muzzin, A., van Dokkuin, P., Kriek, M., Labbe, I., Cury, I., 

Marchesini, D., & Franx, M. 2010, ApJ, 725, 7 
Narayanan, D., Cox, T. J., Hayward, C, & Hernquist, L. 2010, 

arXiv: 1005.3020 
Noeske, K. G., et al. 2007, ApJ, 660, L47 
Nordon, R., et al. 2010, A&A, 518, L24 
Pascale, E., et al. 2008, ApJ, 681, 400 
Pope, A., et al. 2006, MNRAS, 370, 1185 
Perez-Gonzalez, P. G., et al. 2005, ApJ, 630, 82 
Reddy, N. A., Erb, D. K., Pettini, M., Steidel, C. C., & Shapley, 

A. E. 2010, ApJ, 712, 1070 
Rex, M., et al. 2010, A&A, 518, L13 
Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., Perez-Gonzalez, 

P. G., Blaylock, M., Donley, J. L., & Marcillac, D. 2009, ApJ, 

692, 556 
Rigby, J. R., et al. 2008, ApJ, 675, 262 
Rigopoulou, D., Spoon, H. W. W., Genzel, R., Lutz, D., 

Moorwood, A. F. M., & Tran, Q. D. 1999, AJ, 118, 2625 
Roche, P. F., Aitken, D. K., Smith, C. H., & Ward, M. J. 1991, 

MNRAS, 248, 606 
Roussel, H., Sauvage, M., Vigroux, L., & Bosma, A. 2001, A&A, 

372, 427 
Royer, P., et al. 2010, A&A, 518, L145 
Rujopakarn, W., Rieke, G. H., Eisenstcin, D. J., & Juneau, S. 

2011, ApJ, 726, 93 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 
Sanders, D. B., Mazzarella, J. M., Kim, D.-C., Surace, J. A., & 

Soifer, B. T. 2003, AJ, 126, 1607 
Sanders, D. B., et al. 2007, ApJS, 172, 86 
Shi, Y., Rieke, G., Lotz, J., & Perez-Gonzalez, P. G. 2009, ApJ, 

697, 1764 



Siebenmorgen, R., & Kriigel, E. 2007, A&A, 461, 445 

Smith, J. D. T., et al. 2007, ApJ, 656, 770 

Soifer, B. T., Sanders, D. B., Madore, B. F., Neugebauer, G., 

Danielson, G. E., Elias, J. H., Lonsdale, C. J., & Rice, W. L. 

1987, ApJ, 320, 238 
Swinbank, M., et al. 2010a, arXiv:1002.2518 
Swinbank, A. M., et al. 2010, Nature, 464, 733 
Symeonidis, M., Page, M. J., Seymour, N., Dwelly, T., Coppin, 

K., McHardy, I., Rieke, G. H., & Huynh, M. 2009, MNRAS, 

397, 1728 
Takagi, T., et al. 2010, A&A, 514, A5 
Tacconi, L. J., et al. 2010, Nature, 463, 781 
Targett, T. A., Dunlop, J. S., McLure, R. J., Best, P. N ., 

Cirasuolo, M., & Almaini, O. 2010, arXiv:1005.5176 
Thompson, T. A., Quataert, E., & Murray, N. 2005, ApJ, 630, 167 
Thompson, T. A. 2009, Astronomical Society of the Pacific 

Gonference Series, 408, 128 
Tielens, A. G. G. M. 2008, ARA&A, 46, 289 
Truch, M. D. P., et al. 2009, ApJ, 707, 1723 

Veilleux, S., Kim, D.-G., & Sanders, D. B. 2002, ApJS, 143, 315 
Wilson, G. W., et al. 2008, MNRAS, 390, 1061 
Wu, H., Gao, C., Hao, C.-N., Liu, F.-S., Wang, J.-L., Xia, X.-Y., 

Deng, Z.-G., & Young, C. K.-S. 2005, ApJ, 632, L79 
Younger, J. D., et al. 2008, ApJ, 688, 59 
Younger, J. D., et al. 2010, MNRAS, 984 
Y uma, S., Ohta, K Yabe, K., Kajisawa, M., & Ichikawa, T. 2011, 

larXiv:1105.2598l 



